Improved quantum algorithms for the ordered 
search problem via semidefinite programming 
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One of the most basic computational problems is the task of finding a desired item in an ordered 
list of N items. While the best classical algorithm for this problem uses logj A'^ queries to the list, a 
quantum computer can solve the problem using a constant factor fewer queries. However, the precise 
value of this constant is unknown. By characterizing a class of quantum query algorithms for ordered 
search in terms of a semidefinite program, we find new quantum algorithms for small instances of the 
ordered search problem. Extending these algorithms to arbitrarily large instances using recursion, 
we show that there is an exact quantum ordered search algorithm using 41oggQ5 A'' « 0.433 logjAf 
queries, which improves upon the previously best known exact algorithm. 



I. INTRODUCTION 

The ordered search problem (OSP) is the problem of 
finding the first occurrence of a target item in an ordered 
list of items subject to the promise that the target item 
is somewhere in the list. Equivalently, we can remove the 
promise by viewing the OSP as the problem of finding the 
earliest insertion point for a target item in a sorted list 
ol N — 1 items. The OSP is ubiquitous in computation, 
not only in its own right, but also as a subroutine in 
algorithms for related problems, such as sorting. 

To characterize the computational difficulty of the or- 
dered search problem, we are interested in knowing how 
many times the list must be queried to find the location 
of the target item. The minimal number of queries re- 
quired to solve the problem in the worst case is known 
as its query complexity. Using information theoretic ar- 
guments, one can prove that any deterministic classical 
algorithm for the OSP requires \\0g2 ^1 queries. This 
lower bound is achieved by the well-known binary search 
algorithm 1]. 

Quantum computers can solve the ordered search prob- 
lem using a number of queries that is smaller by a con- 
stant factor than the number of queries used in the binary 
search algorithm. The best known lower bound, proved 
by H0yer, Neerbek, and Shi, shows that any quantum 
algorithm for the OSP that is exact [i.e., succeeds with 
unit probability after a fixed number of queries) requires 
at least (IniV - l)/7r « 0.221 log2 N queries 0- In other 
words, at most a constant factor speedup is possible. 
The best published exact quantum OSP algorithm, ob- 
tained by Farhi, Goldstone, Gutmann, and Sipser, uses 
3[log52Af] ~ 0.526 log2iV queries, showing that a con- 
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slant factor speedup is indeed possible |^]. However, 
there remains a gap between the constants in these lower 
and upper bounds. Since the OSP is such a basic prob- 
lem, it is desirable to establish the precise value of the 
constant factor speedup for the best possible quantum 
algorithm: this constant is a fundamental piece of in- 
formation about the computational power of quantum 
mechanics. 

In this article, we study the query complexity of the 
ordered search problem by exploiting a connection be- 
tween quantum query problems and convex optimization. 
Specifically, we show that the existence of an algorithm 
for the OSP that is translation invariant (in the sense 
of jj]) is equivalent to the existence of a solution for 
a certain semidefinite program (SDP). By solving this 
semidefinite program numerically, we show that there is 
an exact quantum query algorithm to search a list of size 
N = 605 using 4 queries. 

Since the size of the semidefinite program increases 
as we increase N , we cannot directly perform a numer- 
ical search for a quantum ordered search algorithm for 
arbitrarily large problem instances. However, by ap- 
plying the 4-query algorithm recursively, we see that 
there is an exact algorithm for a list of size N using 
41oggQ5 N K, 0.433 log2 N queries. Thus, our result nar- 
rows the gap between the best known algorithm and the 
lower bound of 0. In particular, this shows that the 
quantum query complexity of the OSP is strictly less than 
log2 Vn , which one might have naively guessed was the 
query complexity of ordered search by analogy with the 
unordered search problem, whose quantum query com- 
plexity is e(yiv)lll3. 

In addition to providing a way of searching for algo- 
rithms, the semidefinite programming approach has the 
advantage that a solution to the dual SDP provides a cer- 
tificate of the non-existence of an algorithm. Thus we are 
able to provide some evidence (although not a proof) that 
N = 605 is the largest size of a list that can be searched 
with k = A queries, by showing that no algorithm exists 
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for iV = 606. 

The remainder of the article is organized as follows. 
In Section |nl we describe the class of translation invari- 
ant algorithms that we focus on and summarize known 
results about such algorithms. In Section ITTll we show 
how these algorithms can be characterized as the solu- 
tions of a semidefinite program. Finally, in Section IIVI 
we present the results obtained by solving this semidefi- 
nite program, and conclude with a brief discussion. 



II. TRANSLATION INVARIANT QUANTUM 
ALGORITHMS FOR ORDERED SEARCH 

A. Query models for ordered search 

In the standard query model for the ordered search 
problem, sometimes known as the comparison model, a 
query to the xth position of the list indicates whether the 
target item occurs before or after (or at) that position. If 
the target item is at position j e {0, 1, . . . , iV — 1}, then 
its location is encoded in the function fj : {0, 1, . . . ,N — 
1} {±1} defined as 




(1) 



When searching an explicit list with no information 
about its structure other than the fact that it is ordered, 
this function captures essentially all the information that 
is available from examining a given position in the list. 
The ordered search problem is to determine j using as 
few queries to fj as possible. 

The ordered search problem has a kind of symmetry: 
for j g {0, 1, . . . , — 2}, if we change the target item 
from j to J + 1, then we find 



-1 



x = 

1< X < N. 



(2) 



Unfortunately, we must handle what happens at the 
boundary (namely, at x = 0) as a special case. How- 
ever, as observed in 3], we can remedy the situation by 
extending fj to the function gj : 'L/2N ~> {±1} defined 
as 



for {0,1,.. 



\fj{x) 0<x<N 
~\~l3{x-N) N<x<2N 

,N -I}, and 



(3) 



(4) 



for j e {N,N + 1, . . . , 2N - 1}, where all arithmetic is 
done in Z/2N, i.e., modulo 2N. The advantage of using 
this modified function is that the symmetry expressed in 
(0) now appears without special boundary conditions as 
a translation equivariance in the group Z/2N, namely as 



Although the functions gj are defined for all j G Z/2iV, 
it is sufficient to consider the problem of determining j 
with the promise that j e {0, 1, ... , A^— 1}. (Indeed, with 
the quantum phase oracle for gj that we define in the 
next section, it will turn out that j is indistinguishable 
from j + N.) For this problem, the functions fj and gj 
are equivalent, in the sense that any algorithm using one 
type of query can be mapped onto an algorithm using 
the same number of the other type of query. A single 
query to fj can be simulated by simply querying gj on 
the original value oi x G {0, 1, . . . , iV — 1}. On the other 
hand, one query to gj can be simulated by a query to 
fj pre- and post-processed according to Thus, there 
is no loss of generality in using the modified function gj 
instead of the original function fj : the query complexity 
of the OSP is the same using either type of query. 



B. Quantum query algorithms 

In the quantum mechanical version of the query model, 
access to the query function is provided by a unitary 
transformation. Specifically, we will use the phase oracle 
for gj , a linear operator Gj defined by the following action 
on the computational basis states {\x) : x e Z/2iV}: 



gjix)\x) 



(6) 



A fc-query quantum algorithm is specified by an initial 
quantum state \ipo) and a sequence of (j-independent) 
unitary operators Ui,U2, ■ ■ ■ ,Uk. The algorithm begins 
with the quantum computer in the state |V'o)i and query 
transformations and the operations Uj are applied alter- 
nately, giving the final quantum state 



,) ■.= UkG,Uk-i...UiGj\yjo). 



(7) 



We say the algorithm is exact if {4>j\4>j') = 6jj' for all 
e {0, 1, . . . , iV — 1}, since in this case there is some 
measurement that can determine the result j mod N 
with certainty. (Note that the final unitary Uk has no 
effect on whether the algorithm is exact, but it is conve- 
nient to include for the purposes of the following discus- 
sion.) For each value of N, our goal is to find choices of 
\ipo) and Ui,U2, . . . ,Uk for k as small as possible so that 
the resulting quantum algorithm is exact. 

The search for a good quantum algorithm for the OSP 
can be considerably simplified by exploiting the transla- 
tion equivariance lO of the function gj 3]. This equiv- 
ariance manifests itself as a symmetry of the query op- 
erators. In terms of the translation operator T defined 
by 

T\x):=\x + 1) \/x€l/2N (8) 
(where addition is again performed in Z/27V), we have 



gj+i{x) = gj{x - e) yj, X, £ G Z/2N. 



(5) 



TGjT-^ = Gj+i V j e Z/27V. 



(9) 
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Thus, it is natural to choose the quantum algorithm to 
have the translation invariant initial state 



iV'o) 



1 



2N-1 



'2N 



(10) 



x=0 



satisfying T I -00) = 
operations Ut, i.e., 



|?/'o)i and translation invariant unitary 
unitary operators satisfying 



TUtT- 



Ut 



(11) 



for i e {1, 2, . . . , fc}. Of course, while © holds for aU 
j e Z/2N, we are promised that j G {0, 1, . . . , TV — 1}. 
Correspondingly, we can require the N possible orthogo- 
nal final states to label the location of the marked item 
as follows: 



\J 
\J 



N)) 
N)) 



k even 
k odd 



(12) 



(where the separation into k even and odd is done for rea- 
sons explained in 3]). Overall, we refer to an algorithm 
with the initial state I|1U|I . unitary operations satisfying 
and the final states 1)12(1 as an exact, translation 
invariant algorithm (in the sense of Q)- 

An advantage of a translation invariant algorithm for 
ordered search is that, if it can find the target item when 
j — 0, then it can find the target item for all values of j. 
Using ©, we have T-iGjT^ = Gq. Thus 

10,) - {T^UkT-^)Gj . . . {PUiT-^)G,{T^i^o)) (13) 

= TWk{T-'G,T')Uk^i . . . Ui{T-'G,T')\4'o) (14) 

= TWkGoUk-i...UiGo\iko) (15) 

= T=\M- (16) 

In other words, if we find an algorithm whose final state 
in the case j = is given by H12|l , then the final state will 
also be given by (O for j e {1, 2, . . . , TV - 1}. 

C. Characterizing algorithms by polynomials 

Another advantage of translation invariant quantum 
algorithms for the OSP is that they have a convenient 
characterization in terms of Laurent polynomials. A Lau- 
rent polynomial is a function Q : C — *■ C that can be 
written as 



D 



Q{z) - 



(17) 



for some positive integer D, where each qi E C. We 
call D the degree of Q{z). We say Q{z) is nonnegative 
if, on the unit circle \z\ = 1, Q{z) is real- valued and 
satisfies Q{z) > 0. Note that for \z\ — 1, z* = z^^, 
so Q{z) is real- valued on the unit circle if and only if 
q, = gl, for aU i G {0,1,..., D}. If Q(z) = Q{z-^) 



for all z € C, i.e., if qi = q-i for all i € {1, 2, . . . , D}, 
we say Q{z) is symmetric. Thus, Q{z) is nonnegative 
and symmetric if and only if qi = q^i G M for all z G 
{0,1,. ..,£)}. An example of a nonnegative, symmetric 
Laurent polynomial that is relevant to the ordered search 
problem is the Hermite kernel of degree A'^ — 1, 



N-l 



N 



1 

N 



- 1 



z - 1 



(18) 
(19) 



The following result of Farhi, Goldstone, Gutmann, 
and Sipser characterizes exact translation invariant algo- 
rithms for the ordered search problem in terms of Laurent 
polynomials. 

Theorem 1 (01 )■ There exists an exact, translation in- 
variant, k-query quantum algorithm for the N -element 
OSP if and only if there exist nonnegative, symmetric 
Laurent polynomials Qo{z), . . . ,Qk{z) of degree N — 1 
such that 

Qo{z)=Hn{z) (20) 

Qt{z) = Qt-i{z) atz^ = {-lY 

\fte{l,2,...,k} (21) 
Qkiz) = 1 (22) 

— / Qtie'^)dij = l Vt G {0,l,...,fc}. (23) 



Each polynomial Qt{z) in this theorem represents the 
quantum state of the algorithm after t queries. Indeed, 
if we write 



(24) 



i=-(A'-l) 



then 



(zf) =2^(Vt|TV-m)(7V-m-z|0't), (25) 



m— 1 



where 



IV-t) :-[/tGo[/t-i...C/iGo|^o) 



(26) 



is the state of the quantum computer after t queries when 
the target item is j = |3] ■ Given polynomials satisfying 
(|20H23|) . one can reconstruct all of the unitary operators 
Ut for the algorithm using 

Figure shows the (unique) solution to for 
k ~ 2 and A^ = 6 0. In general, the polynomial Qo{z) 
(the Hermite kernel) characterizes complete ignorance of 
the target location at the beginning of the algorithm, and 
subsequent polynomials become flatter and flatter until 
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FIG. 1: Qtie'") as a function of 6* for fc = 2 and iV = 6. The 
solid, long dashed, and short dashed lines represent t = 0, 1, 
and 2, respectively. The intersections at roots of 1 and —1 
are indicated by circles and squares, respectively. 



the final polynomial Qk{z) = 1 is reached, corresponding 
to exact knowledge of the target location. Because each 
query can only change the quantum state in a restricted 
way, successive polynomials must agree at certain roots 
of ±1. Also, each polynomial must be nonnegative and 
suitably normalized. 

With fc = 2, there is a unique choice for the polynomial 
Qi{z), which might or might not be nonnegative depend- 
ing upon the value of TV. For N < 6, this polynomial is 
nonnegative (showing that an ordered list of size < 6 
can be searched in two quantum queries), whereas for 
TV > 7, it is not 0. 

The best ordered search algorithm discovered by Farhi 
et al. was found by considering fc = 3 queries. For fixed 
values of the degree iV — 1, they numerically searched for 
polynomials Qi{z),Q2{z) satisfying the constraints (PUl - 
I23|l of Theorem n The largest value of N for which they 
found a solution was N = 52. Applying this 52-item or- 
dered search algorithm recursively gives an algorithm for 
instances with A^ arbitrarily large. Specifically, one di- 
vides the list into 52 sublists and applies the algorithm to 
the largest (rightmost) item of each sublist, finding the 
sublist that contains the target in 3 queries. This pro- 
cess repeats, with every 3 queries dividing the problem 
size by 52, leading to a query complexity of 3[log52 A^]. 
(Note that although the base algorithm in this recursion 
is translation invariant, the scalable algorithm generated 
in this way is not.) 

In general, recursion can be used to turn small base 
cases into scalable algorithms, so improved quantum al- 
gorithms for the OSP can be found by discovering im- 
proved base cases. Subsequent work by one of us (AJL) 
and collaborators sought such algorithms using a conju- 
gate gradient descent search for the polynomials Qt{z) 
[g. This method is guaranteed to work (for a small 



enough step size) because the space of polynomials sat- 
isfying H2UI - I23|I is convex. The best solutions found by 
this method were A^ = 56 for fc = 3 and A^ ~ 550 for 
k — A, implying a 41og55QA^ « 0.439 logjA^ query re- 
cursive algorithm. Unfortunately, conjugate gradient de- 
scent (or any approach based on local optimization) can 
never prove that finite instance algorithms do not exist 
for a given number of queries k. It could always be the 
case that lack of progress by a solver is indicative of in- 
adequacies of the solver {e.g., the step size is too large, 
etc.). In the next section, we recharacterize exact trans- 
lation invariant quantum OSP algorithms in a way that 
allows either their existence or nonexistence (whichever 
the case may be) to be proved efficiently. 



III. A SEMIDEFINITE PROGRAM FOR 
TRANSLATION INVARIANT QUANTUM 
ALGORITHMS FOR THE OSP 

A. Formulation of the SDP 

In this section, we show that the problem of finding 
Laurent polynomials satisfying the conditions of Theo- 
rem^can be viewed as an instance of a particular kind of 
convex optimization problem, namely a semidefinite pro- 
gram 0. The basic idea is to use the spectral factoriza- 
tion of nonnegative Laurent polynomials to rewrite equa- 
tions H20H23|I as linear constraints on positive semidefi- 
nite matrices. 

The spectral factorization of nonnegative Laurent 
polynomials follows from the Fcjcr-Riesz theorem; 

Theorem 2 Let Q{z) he a Laurent polynomial 

of degree D. Then Q{z) is nonnegative if and only if 

there exists a polynomial P{z) ~ X^iLoP'-^' '^f degree D 
such that Q{z) = P{z)P{l/ z*)* . 

Let Tri denote the trace along the ith super-diagonal 
(or (— i)th sub-diagonal, for i < 0), i.e., for an A^ x Af 
matrix X, 



I 1^1=1 ^e-id i < 0. 



(27) 



The Fejer-Riesz theorem can be used to express nonnega- 
tive Laurent polynomials in terms of positive semidefinite 
matrices, as shown by the following lemma. 

Lemma 1. Let Q{z) = X]iL-^(Af-i) 1^^^ ^ Laurent 
polynomial of degree N — 1. Then Q{z) is nonnegative 
if and only if there exists an N x N Hermitian, positive 
semidefinite matrix Q such that qi = Tr^ Q. 

Proof. The "if" direction follows from the representation 



Q{z) = [1 



-(jv-r)l 



Q 



1 



(28) 
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This Q{z) is real on \z\ = 1 since Q = Q^; it is nonnega- 
tive there because Q is positive semidefinite. 

The converse follows from the spectral factoriza- 
tion of Q{z). Let Q{z) = P{z)P{l/z*)*, let p := 

[pa ■■■ pN-i]'^ , and let z := [l ••• z^"-^]"^. Then 
P{z) — p^z, and Q{z) — z^p*p^z on \z\ — 1. We choose 
Q ■— P*P"^: which by construction is Hermitian and pos- 
itive semidefinite. Furthermore, since Q{z) on \z\ = 1 
determines the coefficients g^, we have Qi = Tr^ Q. □ 

Because the Laurent polynomials in Theorem^are not 
only nonnegative but also symmetric, we can restrict the 
associated matrices to be real symmetric, as the following 
lemma shows. 

Lemma 2. IfQ{z) is a nonnegative, symmetric Laurent 
polynomial, then the matrix Q in Lemma^ can be chosen 
to be real and symmetric without loss of generality. 

Proof. Let Q be a Hermitian, positive semidefinite matrix 
such that Q{z) = z^Qz on \z\ = 1, where z is defined as 
in the proof of Lemma ^ Then the symmetry Q{z) = 
Q{z~^) implies that Q{z) = z^Q^z on \z\ = 1, and by 
averaging these two expressions, we have Q{z) — z^Qz 
on \z\ = 1, where Q := {Q + Q'^)/2 is real and symmetric. 

□ 

Using Lemma 12 we can recast the conditions (|2UII23f) 
of Theorem^as the following semidefinite program: 

Semidefinite Program {S{k,N)). Find real symmet- 
ric positive semidefinite N x N matrices Qq, Qi, . . . , Qk 
satisfying 



Qo = E/N 

% Qt = % Qt-i 

Qk = l/N 
TrQt = 1 



Vie{l,2,...,fc} 
Vi G {0,1,..., fc} 



(29) 
(30) 

(31) 
(32) 



where E is the N x N matrix in which every element is 
1 and 71 : S.^^^ is a linear operator (on the space 

of real symmetric N x N matrices) that computes 
signed traces along the (off-) diagonals, namely 



[Tt X)^ := Tr, X + {-if Tr,_jv X 



(33) 



/orie {1,2,... ,iV-l}. 



The existence of an exact, translation invariant quan- 
tum algorithm for the OSP is equivalent to the existence 
of a solution to this semidefinite program, which can be 
seen as follows: 

Theorem 3. There exists an exact, translation invari- 
ant, k-query quantum algorithm for the N -element OSP 
if and only if S{k, N) has a solution. 

Proof. Given Qo,Qi, ■ ■ ■ ,Qk satisfying S{k,N), let 



0,(2):= [1 ... 



] . Then the 



a nonnegative, symmetric Laurent polynomial; and con- 
ditions (j^ - IS^ imply conditions (pnVE5|) . respectively. 

Conversely, suppose Qo{z),Qi{z), . . . ,Qk{z) are non- 
negative, symmetric Laurent polynomials of degree A^ — 1 
satisfying Let Qo := E/N, let Qk := I/N, 

and let Qi, O2, . . . , Qk~i be positive semidefinite matri- 
ces obtained from Qi{z), Q2{z), ■ . ■ , Qk-i{z) according to 
Lemma H Then (ED and ^ imply ^ and re- 
spectively. □ 

This reformulation of the problem has the advantage 
that semidefinite programs are a well-studied class of con- 
vex optimization problems. In fact, semidefinite pro- 
gramming feasibility problems can be solved (modulo 
some minor technicalities) in polynomial time 0, Il0| . 
Furthermore, there are several widely available software 
packages for solving semidefinite programs 0, IT3. IT3l | . 

Note that by "solving" a semidefinite program, we 
mean not only finding a solution if one exists, but also 
generating an infeasibility certificate (namely, a solution 
to the dual semidefinite program) if one does not. Thus, 
by solving S{k,N) for various values of k and N , not 
only can we extract algorithms from feasible solutions, 
but we can also generate lower bounds for the quantum 
query complexity of the OSP (assuming we restrict our 
attention to exact, translation invariant algorithms). In 
other words, this approach unifies algorithm design and 
lower bound analysis into a single method. 



B. Improved formulation by symmetry reduction 

In moving from the polynomial to the semidefinite pro- 
gramming formulation, we have increased the number 
of real parameters specifying an exact, translation in- 
variant quantum OSP algorithm from {N — l)(fc — 1) to 
N{N -\- l)(fc — l)/2. As benefits, we have put the prob- 
lem in a numerically tractable form, and we are now able 
to prove nonexistence as well as existence of algorithms. 
But the increase in parameters is nevertheless undesir- 
able. 

Fortunately, in our case we can reduce the size of the 
parameter set roughly by half by exploiting symmetry. In 
particular, in terms of the N xN counterdiagonal matrix 
(the counteridentity matrix) 



J 



we have 

Lemma 3. IfQojQi,- 
then so is JQqJ, JQiJ, 







1 

1 



(34) 



,Qk is a solution to S{k,N), 

■ , JQkJ ■ 



symmetry of each matrix Qj implies that each Qj{z) is 



Proof. The matrices JQtJ are positive semidefinite since 
J is unitary. Clearly, JQqJ = Qo and JQkJ = Qk, so 
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(1211 and ijni) are satisfied. Since Tr^ JQtJ = Tr_i Qt by 
the definition of J, and since Tr_i Qt = Tr^ Qt because 
each Qt is a symmetric matrix, is satisfied. Finally, 
(|32|l is satisfied since = I. □ 

Thus, by convexity, if QQ,Qi,---,Qk is a solu- 
tion to S{k,N) then so is \{Q^ + JQqJ),\[Qi + 
JQi J), . . . , \{Qk + JQkJ)- In other words, we can as- 
sume that the matrices Qt commute with J without loss 
of generality. 

We can use group representation theory to harness this 
symmetry. Note that {/, J} is an A^-dimensional (re- 
ducible) representation of the group Z/2. Since J has 
[iV/2j eigenvalues equal to —1 and the rest equal to -1-1, 
this representation can be diagonalized into \N/2\ copies 
of the sign representation and [iV/2] copies of the triv- 
ial representation by some unitary matrix U . The set of 
matrices that commute with all matrices in this represen- 
tation (the representation of the commutant subalgebra of 
Z/2) are therefore block-diagonalizable (by the same ma- 
trix 17) into two blocks, with one block having size [Af/2J 
and the other having size [A^/2] . When N is even. 



U 



1 

71 



J 



(35) 



and when N is odd, 



U 



1 

71 



I I 
V2 
J -J 



(36) 



where / and J arc the [Af/2J by [-/V/2J identity and coun- 
teridentity matrices, respectively. 

Now, since we can choose the matrices Qt to com- 
mute with J without loss of generality, we can block- 
diagonalize them into twice as many matrices, each of 
which has about one quarter the number of elements. 
For example, for N even, Q = JQJ implies that Q has 
the form 



A B 
JBJ JAJ 



(37) 



where A = and B = JB^ J. Thus we have 



1 


I I ' 


T 


A 


B 




7 / " 


2 


J -J 




JBJ 


JAJ 




J -J 


'A + BJ 
















A 


- BJ 







(38) 
(39) 



so that Q is positive semidefinite if and only if ^ ± BJ 
are both positive semidefinite. The net effect of this sym- 
metry reduction is to cut the number of real parame- 
ters in S{k, N) to N{N/2 + l){k - l)/2 (for TV even) or 
{N + lf {k - l)/4 (for N odd), i.e., roughly by half. 



k 


iV* 


2 


6 


3 


56 


4 


605 


5 


> 5000 



TABLE I: Ordered list sizes A''* searchable a fc-query exact, 
translation invariant quantum algorithm such that no such 
algorithm exists for a list of size A'^* + 1. 



IV. RESULTS AND DISCUSSION 

We solved the semidefinite program S{k, N) for var- 
ious values of k and N using the numerical solvers Se- 
DuMi [13, SDPT3 [12, and SDPA [3. These solvers 
use general-purpose primal-dual interior-point methods 
that eventually become limited by machine memory. (Al- 
though there are algorithms for solving SDPs that are not 
based on interior point methods, we did not attempt to 
use such algorithms.) 

The time required to solve S{k, N) was substantially 
reduced by exploiting the symmetry described in Sec- 
tion ^^Bj In addition, it is helpful that the constraints 
are fairly sparse. Nevertheless, we are ultimately lim- 
ited by the fact that the maximum size of a list that can 
be searched increases exponentially with the number of 
queries, so that we can only consider fairly small values 
oik. 

For each fc < 4, we found the smallest value N* such 
that S{k,N*) has a solution but S{k,N* + 1) does not. 
Although we were able to find solutions to 5(5, N) for 
some values of N, we ran out of machine memory before 
we could find an infeasibility certificate. A summary of 
the values N* we obtained is presented in Tabled 

By recursion, the k = A, N* = 605 query algorithm 
yields a scalable algorithm whose query complexity is 



41oggo5 AT « 0.433 log2 A^. 



(40) 



This result also implies improvements to other algo- 
rithms; for example, it implies a quantum sorting algo- 
rithm whose query complexity is AN loggQ5 N. 

As mentioned in the introduction, infeasibility of 
S{k, N* + 1) does not necessarily imply that A^* is the 
largest size of a list that can be searched with a fc-query 
exact, translation invariant algorithm. However, it seems 
reasonable to conjecture that this might be the case. In- 
deed, for fc = 2 and 3, we know that the values of A^* 
in Table are optimal. For those smaller problems, we 
were able to numerically solve the larger SDP developed 
by Barnum, Saks, and Szegedy to characterize general 
quantum query algorithms |l4| . (To solve large enough 
instances of that SDP, it was also crucial for us to exploit 
the symmetry of the problem.) Those results show that 
A^ = 6 and A^ = 56 are the largest sizes of lists that can 
be searched with fc = 2 and fc = 3 queries, respectively. 
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FIG. 2: Laurent polynomial coefficients g,'*' as a function of 
i for exact, translation invariant OSP algorithms. From top 
to bottom, fe = 2, 3, and 4, with A'^ = 6, 56, and 605, respec- 
tively. 



even when the assumption of translation invariance is re- 
moved. 

Whether the In lower bound on the query com- 
plexity of the OSP can be saturated remains open. How- 



ever, the structure of the algorithms we obtained suggests 
the possibility of a well-behaved analytic solution, and it 
would be interesting to understand the behavior of the 
solution in the limit of large N . Figure |21 shows the co- 
efficients of the polynomials Qt{z) associated with the 
optimal feasible solutions Qt for k — 2,3,4. Note the 
similarity of the coefficients for different values of N. 

Not only have we found a particular quantum algo- 
rithm for the ordered search problem, but we have also 
demonstrated the usefulness of semidefinite programming 
as a numerical technique for discovering quantum query 
algorithms. Indeed, the connection between quantum 
query complexity and convex optimization is not unique 
to the ordered search problem: as mentioned above, ar- 
bitrary quantum query problems can be characterized in 
terms of semidefinite programs |14| . Thus, semidefinite 
programming appears to be a powerful tool for studying 
quantum query complexity. 

After this work was completed, we learned that Ben- 
Or and Hassidim have developed a new approach to 
quantum algorithms for ordered search based on adap- 
tive learning p^ . Their resulting algorithm is not exact, 
but rather, is zero-error, with a stochastic running time 
(sometimes referred to as a Las Vegas algorithm). The 
expected running time of their algorithm is 0.32 log2 N. 
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